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A land cover map from Ellesmere Island to Tierra del Fuego: the 
combination of two continental land cover mapping projects and 
comparison to global maps 
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ABSTRACT 


Land cover affects and is affected by climate 
change and therefore is a critical input to climate, 
biogeochemical,  hydrological, and species 
distribution models. There are multiple global 
land cover datasets such as GLC2000, the 
MODIS land cover product, and GlobCover. 
Although all are useful for modeling at the global 
scale, the spatial resolution, thematic detail, and 
accuracy 1s often insufficient for regional studies. 
The North American Land Change Monitoring 
System (NALCMS) was founded in 2006 with 
the goal to provide annually updated land cover 
information for the North American continent, 
which 1s Canada, the United States, and Mexico. 
As a first effort a continental map of the land 
cover in 2005 based on 250m MODIS data was 
completed. The Latin American network to 
monitor and study natural resources (SERENA) 
was founded in 2008 with the aim to monitor, 
study, and distribute information associated with 
changes due to burnt biomass as well as land 
cover of Latin America and the Caribbean. The 
network is formed by 32 participants from 11 
countries and 18 institutions distributed over the 
whole region. The recently completed 2008 land 
cover map is based on 500m MODIS data. This 
study combines both products and presents a land 
cover map from Ellesmere Island to Tierra del 
Fuego. Critical steps in map combination 
included legend  harmonization, adequate 


selection of projection and resolution parameters, 
and adjustments for differences among the years. 
The new map was compared to an existing land 
cover product, the MODIS land cover map of the 
year 2005. For important classes in terms of their 
total area the spatially-explicit correspondence 
among the maps is moderate (approximately 
70%). Other classes such as mixed forest and 
erassland indicate poor agreements due to their 
intrinsic heterogenelty. These dissimilarities may 
be attributed to the methods (supervised or 
hybrid) and data but also to remaining differences 
in legends and processing issues to enable map 
comparison. 


Keywords: ¡Land cover, Map comparison, 
MCD12Q1, NALCMS, SERENA,  MODIS, 
America 


INTRODUCTION 


Land cover and land use maps are probably 
the most common  derivatives of satellite 
images. Land cover maps are required by the 
modeling community for climate, biodiversity, 
and partly also hydrology and landscape analysis. 
Another field of people interested in land cover 
maps the global change community. Land cover 
1s defined at the physical material at the Earth 
surface. Typical materials are grass, trees, barren 
ground, or water. In this sense 1t is different from 
land use that 1s the description of the human use 


of land. Most maps are in reality a mix of both 
and thus termed a land use and land cover map 
(Fisher et al. 2005). Prominent examples are 
the mapping of land use classes urban area 
and agricultural land while other classes 
describe typical land covers, e.g. grassland and 
water. Land cover maps can be derived from 
satellite images in a fairly automatic manner and 
general land uses as described above may be 
inferred with reasonable accuracy. 

Maps can be derived from various spatial 
resolutions, map extents and thematic detail that 1s 
the number of classes. Usually there is 
proportional relationship with increasing spatial 
resolution, decreasing extent and more detailed 
thematic classes. This 1s directly liked to 
feasibility. At least until nowadays there 1s no 
global map at 30m resolution but there exist 
multiple maps between 300 and 1000m. On the 
other hand 1t makes little sense to map a small city 
at that coarse resolution. The thematic detail 
not exclusively refers to the number of classes 
but also to the detail that the classes characterize. 
This paper will focus on the combination of two 
existing land cover — and in the strict sense used 
above — also land use maps for the American 
continent. The maps originate from two existing 
mapping projects in which the authors 
participate. The North American Land Change 
Monitoring System (NALCMS) is a tri-national 
initiative Of Canada (Canada Centre for Remote 
Sensing / Natural Resources Canada — 
CCRS/NRCan), the United States of America 
(United States Geological Survey — USGS), and 
Mexico (National Institute for Statistics and 
Geography — INEGI, National Commission for 
the Knowledge and Use of Biodiversity — 
CONABIO, National Forestry Commission — 
CONAFOR). The work 1s supported by the 
Commission for Environmental Cooperation 
(CEC) that is a tri-national institution to support 
environmental studies for the continent. The aims 
of NALCMS are to develop a land cover 
monitoring system at the continental scale using 
satellite images. It mainly employs images of 
the Moderate Resolution Imaging 
Spectroradiometer (MODIS) sensor. Using 
algorithms developed at CCRS seven reflective 
bands important for land cover monitoring are 
filtered for clouds and errors, downscaled to 
250m spatial resolution and aggregated to 
monthly composites (Latifovic et al, 2012). In 
2010 a baseline map, the land cover of North 
America for the year 2005, was completed and 
made publically available as wall map and in 
digital format (CEC 2010). The map has 19 
thematic classes that characterize different types 
of forest, shrubland, grassland, lichen and 
moss, wetland, cropland, barren land, urban 
area, water, and snow / ice. The overall accuracy 
of the continental map 1s 68% for not 
considering ambiguity in the labeling of reference 


data and 82% 1f an alternative label is considered 
as correctly classified (Latifovic et al, 2012). 
Current activities of this on-going mapping 
initiative are the detection of change area and 
the characterization of changes for map 
updating. The long-term goal is an annually- 
generated map product. However, also 
continuous products and higher spatial resolution 
maps are anticipated. 

Within the NALCMS initiative there are 
common guidelines for product development 
such as the baseline dataset, the legend that 
follows the Land Cover Classification System 
(LCCS, di Gregorio 2005) standards, and partly 
the algorithm. However, each country develops 
1ts own dataset also considering the national 
needs by not losing the  overarching 
international requirements. Specific national 
requirements can be a thematically more detailed 
map that 1s aggregated for the continental map or 
a particular mapping approach to represent the 
diversity of the land cover. The latter was 
necessary for land cover mapping in Mexico 
where in many cases exists more than one label 
for a pixel of 6.25ha (Colditz et al. 2012). Using 
a multitude of discrete classifiers, C5.0 decision 
trees (Quinlan et al. 1993) the proportion of each 
land cover class was estimated. These class 
memberships represent  commonly known 
ecological transition zones (ecotones) due to 
elevation changes or scarcity of water and 
heterogeneous areas with many mixed pixels, 
because the patchiness of the actual land cover 
1s well beyond the resolution of 250m pixels 
(Colditz et al. 2011, Colditz et al 2012). 
Membership maps can be transformed to discrete 
land cover maps by assigning to each pixel the 
class with the highest proportion. 

The same relatively new but computationally 
expensive approach was employed in the project 
“Latin American network to monitor and 
study natural resources” (SERENA). Founded 
in 2008 and embedded in the regional Global 
Observation of Forest and Land Cover Dynamics 
(GOFC-GOLD) program  RedLaTIF, the 
network brings together 52 researchers from 10 
countries distributed over the whole Latin 
American continent. The aim of this initiative 1s 
to monitor, study, and distribute information 
associated with changes due to burnt biomass and 
land cover from Mexico to Tierra del Fuego. The 
analysis 1s based on a set of monthly composites 
of MODIS images with a spatial resolution of 
500m filtered for clouds and processing issues. 
The land cover sub-group recently finished the 
land cover map of Latin America for 2008 
(Blanco et al. 2012). The wall map was released 
in spring 2012, and the digital release 1s 
anticipated for the end of 2012. The overall map 
accuracy was 84% on the continental scale and 
major confusions in the error matrix are related 
to ecological gradients and less to climate zones. 


Both, NALCMS and SERENA maps are 
comparable because a similar legend was used in 
both projects. The map of 19 classes in the 
continental NALCMS was amplified to 22 classes 
in SERENA. The main difference is explicitly 
mapping sub-tropical classes in SERENA. In 
NALCMS there are combined tropical and sub- 
tropical classes, because their distinction was not 
imperative due to the small area proportion of 
their extent. In SERENA, however the majority 
of the mapping area is located within the tropical 
and sub-tropical climate zones. 

The following chapter will describe how the 
maps were combined and which adjustments 
were necessary to link classes thematically. The 
results chapter will present the final map for the 
American continent and will compare this map 
quantitatively to an existing global land cover 
map. The final chapter outlines the conclusions. 


2. METHODS FOR MAP 
COMBINATION AND ANALYSIS 


The NALCMS as well as the SERENA maps 
were generated in the Lambert Azimuthal Equal 
Area (LAEA) projection but with two different 
projection origins. Also, the NALCMS map has a 
spatial resolution of 250m in contrast to 
SERENA with 500m. For the generation of the 
America map both maps were reprojected to a 
new projection center (1095, 902W) and a spatial 
resolution of 500m. 

A second step concerned legend merge of the 
19 classes NALCMS and 22 classes in 
SERENA. 

The main difference between legends of both 
maps concerns the sub-tropical classes that 
were discerned in SERENA but not in NALCMS. 
However, the issue can be solved, because 
tropical classes in the NALCMS map extent 
only exist im Mexico, and Mexico was also 
covered by SERENA. Therefore the new 
America map employed data of SERENA for 
Mexico. 


NALCMS data were only used for Canada and the 
United States where all tropical and sub-tropical 
classes were labeled “sub-tropical”. Other classes 
of lichens and mosses in NALCMS only exist 
in northern Canada and Alaska and not at all 
in South America. 

Quantitative comparison of the America map to 
the MCD12Q1 map to another existing map was 
based on the comparison matrix. First, both maps 
were transformed to a common generalized 
legend of basic land cover classes. Cells of the 
diagonal of the matrix indicate the agreement (A) 
between both maps (Equation 1). In contrast, 
disagreement (D) can be calculated for both maps 
(X, Y) and 1s the sum of the column or row totals 
(+k, k+) minus the diagonal cell (Equations 2 
and 3). An additional division by 2 is 
necessary for not double counting commission 
and omission among different classes. 


A(X Y), = N kk (1) 


D(XX)h _"+ Kk — Try 
2 (2) 


DY): ¿E TEE 
2 (3) 


3. RESULTS AND ANALYSIS 


This section will present the results of this study. 
First, the land cover map for the entire American 
continent from Ellesmere Island to Tierra del 
Fuego will be presented and described. 
Another section will focus on a quantitative 
analysis of this map and a comparison to another 
existing map product for the American continent. 
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3.1. The America map 
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Figure 1. Land cover map of the American continent. 


Figure 1. Land cover map of the American continent 





Figure 1 presents the land cover map for the 
American continent. Zonal patterns dominate the 
northern parts of the North American continent. 
Meridianal mountainous systems like the 
Appalachianss and the Western Cordilleras 
interrupt this pattern with broadleaf deciduous 
forests and shrubland in the semi-arid regions 
of the interior basins and Colorado plateau. 
Also, land use patterns become visible at a broad 
scale with cropland extending from Alberta and 
Saskatchewan to the southern United States 
East of the Mississippi/Missouri system. West 
of this system dominate temperate grasslands. 
Towards Mexico those classes change towards 
their  tropical-subtropical counterparts with 
temperate classes continuing towards the South at 
higher elevations. Other Central American 
countries are dominated by  cropland that 
intermingles with tropical and  sub-tropical 
forests. South America is dominated by the 
tropical broadleaf evergreen forests in the Amazon 
basin that extent towards the Atlantic coast with 
transitions to tropical shrubland and sub-tropical 
forests. Climatic gradients are clearly visible for 
South  America?s lowlands, changing from 
tropical to sub-tropical and  eventually to 
temperate classes in Patagonia. Also the 
hypsometric gradient can be easily observed for 
land cover classes in the Andean ranges. 
Latifovic et al. (2012) and Blanco et al. (2012) 
describe in more detail and quantitative fashion 
the land cover for North America and Latin 
America, respectively, and Colditz et al. (2012) 
with even more detail for Mexico. 


3.2. Quantitative map comparison 


For quantitative map comparison the North 
American map was aggregated to 12 basic land 
cover classes. The specific aggregation scheme is 
presented in Table 1. This generalized map was 
compared to the MODIS land cover product for 
the year 2005 with 17 classes following the IGBP 
legend (Friedl et al. 2002, Friedl et al. 2010). 


The MODIS product was obtained in tiles that 
were mosaicked for entire America and 
projected to the same LAEA parameters and 
grid as the America map. 

The aggregation rules of the MCD12Q1 product 
with S00m spatial resolution was not stralght- 
forward and followed specific semantic and 
spatial rules that are also summarized in Table 1. 
Figure 2 presents the quantitative class-wise 
comparison between the generalized America 
and MCD12 IGBP map. Coniferous forest that 
mainly exists in Canada, the United States and 
higher elevations for the rest of the study area 
shows a good agreement of 62%. There is a 
proportion of 25% that was only mapped in 
MOD12 IGBP but not in America and a 
proportion of 13% vice versa. The agreement was 
even better for broadleaf forest with 94%. Even 
more important 1s that this class makes up a 
large part of the land surface with approximately 
28% of the total area. In this light one has to see 
the poor agreement of mixed forest (35%) 
because less than 4% of the total land 
corresponds to this diverse class that is hard to 
define and only exists in North and Central 
America. Another spatially important class is 
shrubland with approximately 16% of the total 
land surface. This class shows good 
correspondences above 70%. Grassland indicates 
lower agreements with only slightly above 50%. 
The larger disagreement occurs for the MOD12 
IGBP map with 33% because this map classified 
more grassland than the America map. Lichen 
and moss only exists in the northern zones of 
North America and shows good general 
agreements above 70%. Although having poor 
agreements (23%) this number for wetland has 
only a limited importance because of its small 
area proportion over the entire study area 
(approximately 2%). In contrast, cropland with 
a share between 11 and 15% of the area 1s rather 
important. The agreement of 70% 1s moderate 
due to more mapped cropland in America than in 
the MOD12 IGBP map. Other classes (barren 
land, urban and built-up, water, snow and ice) 
show moderate to poor agreements but their 
importance in terms of total area 1s low. 


Table 1. Ageregation of the America map legend and the MCD12 IGBP legend 
to a generalized 12-class legend. 


Genera! land 

cover classes 

1. Needieleaf | 1, Sub-rropical needleleaf forest, 2, Temperate | 1. Evergreen needieleat forest, 

forest or sub-polar needieleaf forest, 3. Sub-polar 3. Deciduous needieleaf forest, 
taiga needleleaf forest 8. Woody savanna 

2. Broadicaf |4. Tropical broudicat evergreen forest, 5. Sub- |2. Evergreen broadleat forest, 4. 

forest tropical broadleaf evergreen forest, 6, Tropical | Deciduous broadleaf forest, E. 
broadleaf deciduous forest, 7. Sub-tropical Woody sevenna, 9. Savanna, 14. 
broadleaf deciduous forest, E. Temperate Cropland / natural vegetation 
broadleaf evergreen forest, 9. Temperate or mosaico 
sub-polar broadleaf deciduous forest 

3. Mixed 10. Mixed forest 5. Mixed forest 

forest 

4. Shrubland | 11. Tropical shrubland, 12. Sub-tmpical 6. Closed shrublands, 7. Open 
shrubland, 13. Temperate or sub-polar shrublands, $. Woody savanna, 
shrubland 9. Savanna, 14. Cropland / 

natural vegetation mosale 

5. Grassland | 14. Tropical grasstand, 13. Sub-tropical 8. Woody savanna, 9. Savenna, 
grassland, 16. Temperate or sub-polar 10. CGrrasslands, 14. Cropland 
grassland natural vegetation mosaic 

6. Lichen 17. Sub-polar or polar shrubland-lichen-moss, |7. Open shrublands 

and moss 18, Sub-polar or polar grassiland-lichen-M05s, 
19. Sub-polar or polar barren-lichen-mass 

7. Wetdand |20. Wetland 11. Permanent wetlands 

8. Cropland |21. Cropland 12. Croplands, 14. Cropland : 

natural vegetation mosalc 

9. Barren 212. Barren land 7. Open shrublands, 16. Barren 

land or sparsely vegetated 

10, Urban . Urban arca 13, Urban end built-up 

and bullt-up 

11, Water . Water 0. Water 

12. Snow 25. Permanent ice and snow 15. Snow and ice 

and ¡ce 


Ámerica map MCD12 IGBP 





Note: Some classes of the MCD12 IGBP map correspond to multiple land cover classes in the generalized 
legend. In this case spatial rules were applied, 1.e. this class was converted to a specific general land 
cover class 1f 1t corresponded to the label. 


An interesting number can be presented if the other (remaining classes). For forest the 
generalized 12 classes are combined to 3 basic agreement 1s 90%, for natural vegetation 81%, 
classes that are forest (coniferous, broadleaf and and for all other classes combined 76%. For the 
mixed forest), natural vegetation (shrubland, map overall, the agreement is 72%. 


grassland, lichen and moss, and wetland) and 


12. Snow and ice 2 

11. Water | 

10. Urban and built-up | 
9. Barren land 
8. Cropland 

7. Wetland | 


6. Llichen and moss 


5. Grassland 


A. Shrubland 
3. Mixed forest 
2. Broadleaf forest 4 


1. Needleleaf forest 


0% 20% 


A0% 60% 80% 100% 


mM America M America AND GlobCover "MOD12IGBP 





Figure 2. Class-specific agreement and disagreement (in per cent and Mio ha) between the 
America and MDD12 IGBP maps with generalized legends (see Table 1) 


4. CONCLUSIONS 


The paper presented the generation of a land 
cover map for the entire American continent. Due 
to the similarities between the NALMCS and 
SERENA project such a map can be 
accomplished by only a few processing steps. 
The map is important for projection on the 
continental scale such as climate modeling and 
species distribution and migration studies. The 
map 1s similar to other existing products with the 
advantage of presenting more thematic classes 
than several other products. 
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ABSTRACT 


Tropical deforestation is today considered as a 
major environmental issue. And it is accepted 
that deforestation and forest degradation are 
estimated to contribute to about 7 to 153 % of the 
total greenhouse gas emissions. A post-Kyoto 
protocol mechanism was developed, termed the 
Reduction of Emissions from Deforestation and 
forest Degradation (REDD+). The development of 
deforestation monitoring systems and 
determination of historical trends for forest cover 
change are key point in the establishment of 
national or sub- national carbon policy. This paper 
aim to show how through three different projects, 
and four different countries, a single standardized 
methodology has been developed and applied to 
monitor forest cover in 1990, 2000 and 2010. 
Each of those projects (project GSE FM REDD in 
Gabon funded by the European Spatial Agency, 
REDDAF in CAR and Cameroon funded by the 
European Commission and SLCS in Bolivia 
funded by the CNES) aim to make an inventory 
of forest areas at national or sub-national scale in 
order to take stock of the forest cover changes over 
the past 20 years. The proposed services will form 
the basis of a wall-to-wall regional- wide forest 
cover monitoring in this four countries. Therefore, 
a new and cost efficiency methodology was set up 
in order to automatically map forest cover using 
only optical data. 


Keywords: Remote Sensing, Image Processing, 
Forest monitoring. 


1. INTRODUCTION 


Tropical deforestation is today considered as a 
major environmental issue Justice et al. (2001), 
Zhang et al. (2005), Hansen et al. (2008). 
Deforestation and forest  degradation are 
estimated to contribute to about 7 to 15 % of the 
total greenhouse gas emissions Achard et al. 
(2007), GOFC-COLD (011). A post-Kyoto 
protocol mechanism was developed, termed the 
Reduction of Emissions from Deforestation and 
forest Degradation (REDD) and initiated at the 
UNFCCC Conference of Parties (COP) 11 in 
Montreal in 2005. One of the fundamental 
requirements in the UNFCCC process is the 
establishment of robust and transparent national 
forest monitoring systems. Moreover 
determination of historical trends for forest 
cover change is also a fundamental input in the 
establishment of national or sub-national policy in 
relation to REDD. 

Land cover maps fulfill several functions in 
monitoring forest cover change. Firstly, they serve 
as a reference baseline against which future 
change can be assessed. Secondly, they help 
establish forest areas that need to be monitored for 
change. When using a land cover map to assess 
future change, consistent methodology and spatial 
resolution are critical for the interpretation of 
results Achard et al., (2007). Different methods 
are appropriate and reliable for forest cover 
mapping or monitoring at regional or national 
scales using remote sensing data. An approach 
using  spatially explicit land conversion 
information derived from sampling. 


Bodart et al, (2009) or wall-to-wall mapping 
techniques 1s considered the most  appriate 
(GOFC-COLD, 2011). 

Herein, a standardized method for monitoring 
forest cover and forest cover change from optical 


Tablel. Forest cover mapping project undertaken by sirs 


Earth observation data was developed and applied 
over significant areas in four separate countries 
listed in table 1. 


country Funded by Area 


Bolivia (Pando) 
Gabon (1/3 country in 


SLCS Carbone CNES COINACAPA | 63 000 KM? 


2010) 

Central African Republic REDDAF EU Commission MEE RCA 87 000 KM? 
(Sangha Mbaéré, Lobaye 

and Ombella-M”Poko) 


Cameroon 
Province) 


The estimation of deforestation is a function of 
the forest definition, in the absence of national 
forest definition; the service will comply with 
FAO forest definition IPCC (2006), which 
provides a range of values from which 
countries can select their forest definition. 
Forest should be an area of minimum 1 ha, with a 
minimum tree crown cover (or equivalent stocking 
level) of 30% and a minimum tree height at 


(Central | REDDAF EU Commission MEE 90 000 KM? 
CAMEROON 


maturity of 5m. 

The proposed services will form the basis of a 
wall-to-wall regional-wide forest cover 
monitoring. Therefore, a new and cost efficiency 
methodology was set up in order to automatically 
map forest cover using only optical data. 
Technical specifications of forest area maps and 
forest cover change maps for 1990, 2000 and 2010 
are listed in table 2. 


Table 2.Forest Cover Maps 1990, 2000, 2010 and Forest Cover Change Maps 1990-2000 and 
2000-2010 product specifications 


Forest cover forest cover change map 


Criteria Specification Criteria Specification 


Geometric accuracy Around 30 meters Geometric accuracy Around 30 meters 


Reference system WGS 1984 UTM Reference system WGS 1984 UTM 


Image quality 30 metres Image quality 30 metres 
Image acquisition +/- 3 years Image acquisition +/- 3 years 
Minimum mapping 

pe l ha Minimum mapping unit l ha 
Map classes Forest / Non Forest Map classes IPCC 
Thematic accuracy 85% +/- 5% Thematic accuracy 85% +/- 5% 


In Gabon, the ESA GSE FM REDD extension 
project aims at developing a pre- operational 
system for monitoring forest cover within the 
REDD framework. The activities in Gabon are 
coordinated by SIRS and the overall project by 
GAF AG. The Gabonese Agency for Space 
Studies and Observations (AGEOS) was set up in 
2010 with the aim to establish a national 
infrastructure for environmental monitoring and 
preventing the impacts of climate change. One of 


the objectives of the AGEOS is to develop its 
capacity to monitor forest cover at national level. 
The Gabonese authorities see the GSE FM 
REDD project as a precursor to the establishment 
of such national forest monitoring system. 

In Bolivia, the main aim of the SLCS project 
funded by the CNES and coordinated by SIRS is 
to develop a monitoring system to facilitate access 
to payment for environmental services by local 
fair-trade organizations involved in agroforestry 


and deforestation avoidance projects. The case 
study in the Pando district in Bolivia has been 
developed in partnership with an organization 
of producers of Brazil nuts (COINACAPA). 
These nuts are not cultivated but harvested in 
the Amazon rainforest. This case study aims to 
show that it is an activity that contributes to 
preserving forest resources in the area in which 
they operate. 

Finally, the FP7 REDDAF project coordinated 
by GAF AG, seeks to develop and implement 
improved methodologies for assessing activity 
data and emission factors which are key input data 
for carbon accounting in the REDD process. The 
improved methodologies support the development 
of operational service chains for forest 
monitoring. Two  demonstration areas were 
selected to develop and implement these 
methodologies, one in Central African Republic 
(Sangha Mbaéré, Lobaye and Ombella-M”*Poko) 
and a second in Cameroon (Central Province). 
Both service providers SIRS and GAF AG were 
involved in the production of the map products to 
further demonstrate the interoperability of the 
mapping service developed. 


2. METHODOLOGY 


The Forest Cover Map (Forest/Non-Forest map) 
generation is based on Optical imagery and the 
overall production process follows a stepwise 
approach. 


2.1 Data source 

The production of the forest cover and forest 
cover change maps were exclusively based high 
resolution Optical data. For the four study area, in 
1990 and 2000, Landsat data appears to be the 
most suitable data for historical mapping: datasets 
are freely available (USGS), the satellite has a 
repetitivity of 16 days and the spatial resolution 
meets the services “criteria (cf. table 1) Thus, 
Landsat TM 4 and TM 5 were used for 1990 and 
Landsat 7 ETM+ for 2000. For 2010, in Gabon, a 
combination of Landsat 7 ETM+ SLC off and 
ASTER data were used, in Pando (Bolivia) only 
Landsat 7 ETM+ SLC off were used. In CAR 
and Cameroon, a mix of RapidEye and Deimos 
data were used. Due to the heavy cloud cover 
conditions generally experienced over tropical 
area, several images were required to cover a 
scene entirely. It was almost always possible to 
acquire the imagery within the required time 
interval (1.e. +/- 3years from the reference year). 


2.2 Pre-Processing 

All data acquired were already georeferenced by 
the provider to an UTM projection. Theimage-to- 
image registration was performed using Landsat 
GLS-1990/2000/2005 as reference dataset Gutman 
et al. (2008). All input dataset was initially 
checked visually and a minimum set of 10 evenly 
distributed check points were selected to assess the 
quality of the registration. In cases when the 
computed RMS exceeded 30m, the image was 
rectified using a first order polynomial model. All 
were resampled to an output pixel size of 30m to 
create a consistent dataset for change detection. 
Finally a topographic normalization was applied 
for each image. 

A cloud detection algorithm was developed based 
on the combination of an  unsupervised 
classification approach combined with a visual 
comparison of the results with the input image to 
determine the threshold between cloudy and 
non-cloudy pixels. A cloud mask was produced 
for each input image to derive cloud cover 
statistics and to determine the effective coverage. 
Gap stripe on Landsat 7 ETM+ SLC off were 
processed as cloud area and thus were masked and 
considered as no data area. 


2.3 Thematic Processing 

Depending on statistical results of data covers 
calculations for each reference years (1990,2000 
and 2010), the reference year getting the best 
coverage of data is chosen as time O (10). 10 was 
selected as it provided the most complete and 
homogeneous image coverage will be used as the 
basis of classification forest / non-forest. For 
forest cover map and forest cover change map in 
Gabon and Bolivia, 10 correspond to the reference 
year 2000, in CAR and Cameroon, t0 correspond 
to reference year 2010. 

In a first step, the segmentation and classification 
of EO dataset from 10 was carried out, followed 
by manual/visual enhancements of the classified 
Forest Cover Map. The derived classified map 
product from +0 and the Landsat data from t£/ 
are used to derive the Forest Cover Map for tl 
and the Forest Cover Change Map  10/tl 
(classified into IPCC compliant Land Use 
classes). The 10 or t/ Forest Cover Map and the 
current EO data set from 12 are then the basis to 
derive the Forest Cover Map from 12 and the 
respective mapped Forest Cover Changes. A 
more detailed description of the single processing 
step 1s outlined in the following paragraphs. 


2.3.1 Semi-Automatic Classification of Forest 
Cover and Segmentation Process 

The classification of the area of interest into Forest 
$ Non-Forest was performed for each cloud and 
cloud shadow masked image, by using the 
ERDAS Imagine software. Optical data from 10 
were used as input. An unsupervised classification 
of the area of interest was conducted with a large 
number of spectral classes to ensure a good 
representation of the thematic classes” variability. 
An interactive comparison of the classified image 
was performed to label each spectral class 
according to the Forest/ Non Forest thematic 
classes. A post- processing classification routine 
was applied to vectorise the results and eliminate 
polygons smaller than 1 ha. 


2.3.2 Manual Post-Processing 

Subsequent to the automated classification, an 
Intensive interactive post-processing was applied. 
Areas incorrectly classified were  relabelled. 
The manual post-processing was performed 
using in house developed tools built on the 
ESRI software suite, especially designed for this 
type of activity. Additional information sources 
(e.g. topographic maps) can be used to support 
this interpretation. In case of discrepancies 
between the ancillary data sources and the 
satellite images, the most up-to-date information 
from the satellite imagery was used. 


2.3.3 Detection of Forest Cover Changes 

The 10 Forest Cover Map was used as a Mask 
with the t/ imagery to identify areas of change. 
The geometry of the 10 map was kept when no 
change was detected to avoid the creation of 
artefacts. This process was repeated for 12. The 
detected change areas identified as part of the £/ 
and 12 Forest Cover Map production process were 
classified into cropland, grassland,  wetland, 


settlement or other land use which are IPCC 
compliant. 


2.4 Accuracy assessment of classification results 
As shown in Table 1, the thematic accuracy 
should reach 90%. Many methodologies to assess 
the accuracy of a map have been developed 
Stehman, (2009), Czaplewski et al (2004), 
Czaplewski1 (2003), Congalton, (2001), (1999), 
(1991). The main objectives of the 
implementation of the thematic map accuracy 
assessment procedure are to provide objective 
unbiased information on the map accuracy, and a 
cost-efficient and well accepted procedure. The 
easiest way to ensure that a sample is unbiased is 
to select each element randomly Stehman and 
Czaplewski, (1998). Despite the fact that, random 
selection 15 completely unbiased, it can leave 
some areas uncovered resulting in a non-fully 
representative sample. Therefore, a good 
alternative is a combination of systematic and 
random approach resulting in a  randomly 
unaligned systematic sample scheme shown in 
figure 1 Czaplewski et al, (2004), Stehman and 
Czaplewski, (1998). It consist in a systematic grid 
of 20 by 20 km cell within a single sample of 2 by 
2 km is randomly selected. Each segment was 
photo-interpreted by a qualified photo-interpreter 
independent from the production team. Then 
50 points are randomly selected in each sanple. 
A confusion or error matrix is usually produced 
Congalton, (1991) based on a comparison of 
single points on the maps. The overall 
classification accuracy can be calculated together 
with user's and producer?s accuracies which 
represent respectively errors of omission and 
commission. This approach is particularly suited 
for the Forest Cover and Cover Change Map 
products. 














Figure 1. Illustration of the unaligned systematic sampling strategy over Gabon 


3. RESULTS 

Table 2 show initial results of forest cover by will be done during the phase 2 of REDDAF, 
project and the table 3 the associated overall Table 3 resume the  deforestation rate by 
accuracy by reference year. In Cameroon project and confirm the generally low level of 
results must be confirmed by the deforestation in these four regions and in Gabon 
improvements of the dry-forest classification, it and Cameroon in particular. 


Table 2. Forest covers for 1990, 2000 and 2010 


GABON CAMEROON CAR E OLIVIE 


1990 39.88% 13.14% 66.56% 96.40% 
2000 89.64% 72.92% 65.56% 95.80% 


2010 86.13%:* 74,14% 64.03% 94.40% 


* Deforestation rate for Gabon for a third of the country (102 000km?) 


Table 3. Overall accuracy for 1990, 2000 and 2010 


Overall accuracy 1990 98.8% [95.05% 90.1%  P82% 


Table 4. Deforestation rate by country for 1990/2000, 2000/2010 and 1990/2010 


| GABON | CAMEROON CAR... BOLIVIE 


E 26% 
1990/2000 0.209? 0.22% 0 


2000/2010 0.04% -0.52% 1 
1990/2010 0,24% -0.30% 2 


IDeforestation rate for Gabon between 1990 and 2000 are for the whole country (267 000kim?) 
2Deforestation rate for Gabon for a third of the country (102 000kn?) 


The change results between 2000 and 2010 seem to overestimate the increase in forest areas due to 
misclassification of Dry Forest areas in the Centre Province of Cameroon. In phase 2 of REDDAF, the Dry 
Forest thematic will be more investigated to improve the overall results of the Deforestation mapping. 


4. CONCLUSIONS 


It was shown that it was possible to make wall 
to wall forest cover map at three time period at 
national or sub national scale using only optical 
data in tropical area known forheavy cloud 
cover, especially in Gabon. The maps produced 
are fully compliant with the product 
specifications agreed at the beginning of the 
project with all the users. The results shown 
indicate that the methodology developed is 
applicable in a range of different environments. 
The examples of CAR and Cameroon illustrate 
that a homogenous forest cover map can be 
produced by two separate service providers. In 
addition, as 1t 1s based on relatively simple and 
robust methods, it can potentially be implemented 
by countries wanting to develop their own 
forest monitoring system in the framework 
of REDD. And at international level to produce 
standardized reference forest cover map. 
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ABSTRACT 


Schistosomiasis caused by  Schistosoma 
mansoni 1s a disease conditioned to the 
presence of snails of aquatic habits of the 
genus Biomphalaria. Since schistosomiasis is a 
disease determined in space and time by 
environmental variables, the Geographic 
Information System (GIS) and Remote Sensing 
(RS) are techniques for the identifying factors 
and defining areas of possible risk of 
schistosomiasis. The objective of this study 
was evaluate the using of GIS and RS to 
characterize the spatial distribution of 
schistosomiasis in Mosqueiro Island (MD 
focusing on the construction of scenarios 
representing potential areas for the occurrence 
of the disease. B. straminea was the only 
species of Biomphalaria found in 20 breeding 
sites. They are located in eight census tracts, 4 
sectors located in Southwest MI and 4 in the 
West part of the island. The West region has a 
greater precariousness in relation to sanitation 
than the Southwest part of the MI. The results 
using kernel showed that there are 4 clusters of 
Biomphalaria. They are situated in only one 
cluster in Southwest and 3 in West of MI. The 
satellite image showed that the agglomerate of 
the Southwest having a large area with 
possibility of urban expansion. The lack of 
basic sanitation to support the new possible 
urban areas added with the Biomphalaria 
presence can leads to faecal contamination of 


natural aquatic environments, resulting in new 
sources of schistosomiasis transmission. This 
study showed the importance of the use of 
GIS and RS to study the risk of 
schistosomiasis in the MI. 


Palavras-chave: geographic information 
system, schistosomiasis, Biomphalaria, sistema 
de informacOes geográficas, esquistossomose, 
Biomphalaria 


1. INTRODUCAO 


A  esquistossomose  mansónica,  conhecida 
popularmente como barriga d'água é uma 
doenga parasitária de  veiculacáo hídrica 
causada pelo Schistosoma mansoni que se 
aloja nas veias do sistema porta-hepático; tem o 
homem como seu hospedeiro definitivo e 
caramujos de  água doce  (Gastropoda, 
Pulmonata, Planorbidae) do género 
Biomphalaria como hospedeiros intermediários 
(Brasil, 2008). 

É uma doenca considerada de importáncia em 
saúde pública por estar  intimamente 
relacionada ás condicóes de pobreza, fluxo 
migratório, precariedade do saneamento básico, 
dentre outros fatores socioambientais. As 
manifestacóes clínicas da doencga sáo muito 
variáveis e  dependem da fase e dos 
mecanismos  fisiopatogenéticos  envolvidos 
(Bichara et al., 1997). Na fase inicial os 


sintomas mais aparentes decorrem em geral 
da penetracáo de cercarias na pele o que pode 
resultar em um  quadro clínico com 
manifestacóes  pruriginosas  cChamado de 
dermatite cercariana. Já a fase crónica da 
doenca pode-se apresentar várias formas: 
hepatointestinal, hepatoesplénica e neurológica. 
No Brasil a esquistossomose foi introduzida 
com a importacáo de populacóoes africanas em 
regime de escravidáo em meados do século 
XVI para trabalhar nas plantacóes de cana-de- 
acúcar no nordeste brasileiro (Barreto, 1982). O 
país é considerado oO principal foco das 
Américas, com uma estimativa de cerca de 6 
milh0es de indivíduos infectados e outros 25 
milhóes sob o risco de contrair a doenca 
(Katz e Almeida, 2003). As regides n rdeste 
e sudeste do Brasil concentram as áreas com 
endemicidade mais elevada. 

Com o fim do regime de escravidáo, houve um 
expressivo aumento no fluxo  migratório 
interno, o que levou a expansáo dos focos da 
doenca para outras regid0es até entáo indenes. A 
esquistossomose  chegou a  Amazónia 
impulsionada pelo ciclo económico da borracha 
no final do século XIX. O primeiro registro da 
ocorréncia na Amazónia foi feito por Lutz em 
1919, no estado do Acre. Posteriormente, Davis 
(1934), Pará (1948) e outros observaram vários 
casos em Belém e outras localidades, sem, no 
entanto a comprovacáo de sua autoctonia, O 
que só veio a ocorrer de fato em 1951 quando 
da descoberta do primeiro foco de transmissáo 
ativa feita por Machado e Martins (1951), no 
município de Fordlándia, localizado a margem 
direita do rio Tapajós. Posteriormente, foram 
identificados dois outros focos autóctones: um 
no atual município de Quatipuru descrito por 
Pessoa de Melo e Gueiros (1959), que 
permanece ativo até os dias atuais, e o outro no 
bairro do Reduto em Belém, descrito por 
Galváo (1968), tendo se expandido depois para 
outros bairros da capital paraense (Bichara et 
al., 19977). 

De acordo com a lei municipal N* 7682/94, 
os 71 bairros que integram o município de 
Belém estáo distribuídos em oito Distritos 
Administrativos (DA), a saber: DABEL 
(Belém), DABEN (Benguí), DAENT 
(Entroncamento), DAGUA (Guamá), DAICO 
(Icoaraci)) DASAC (Sacramenta), DAOUT 
(Outeiro) e DAMOS (Mosqueiro) (Belém, 
1994), sendo que apenas duas destas áreas 
(DAOUT e DAMOS) sáo indenes para 


esquistossomose, muito embora já tenha sido 
registrada a presenga do caramujo em vários 
bairros (Aeroporto, Chapéu Virado, Ariramba, 
Baia do Sol, Carananduba, Maracajá, Vila, 
Murubira, Paraíso e Sáo Francisco), bem 
como casos importados da doenca já foram 
notificados (Belém, 2008). 

Como a esquistossomose é uma doenca 
determinada no espaco e no tempo por 
variáveis socioambientais, oO Sistema de 
Informacáo Geográfica (SIG) e O 
Sensoriamento Remoto (SR) sáo técnicas para 
identificar esses fatores e definir áreas de 
possível risco de esquistossomose. (Guimaráes 
et al., 2009). 

Neste contexto, faz-se de extrema importáncia 
O levantamento de dados que possam contribuir 
para um «melhor conhecimento acerca da 
dinámica de transmissáo e do comportamento 
epidemiológico do agravo, os quais poderáo 
ser incorporados como práticas de caráter 
preventivo junto ás estratégias desenvolvidas 
pelo Programa de Saúde da Família (PSP). 

O objetivo deste estudo foi avaliar o uso de 
SIG e SR para caracterizar a distribuicáo 
espacial da  esquistossomose no Distrito 
Administrativo de Mosqueiro, Belém/PA com 
foco na construcáo de cenários que representam 
as áreas potencials para a ocorréncia da doenca. 


2. METODOLOGIA DE TRABALHO 


O Distrito Administrativo de  Mosqueiro 
(DAMOS) é uma ilha fluvial situada entre 
01%03” e 01%05” latitude Sul e 4829” e 4818" 
longitude Oeste (W) de Greenwich, distante 
apenas 44,3 km de Belém e ocupando uma área 
de aproximadamente 212 km” de extensáo 
(Figura 1), sendo limitada a sudoeste pela 
Baia do Guajará, a oeste pela Baía de 
Santo António, a noroeste pela Baía do 
Marajó, ao norte e nordeste pela Baía do Sol, ao 
sul pelo Furo do Maguari e a sudeste pelo Furo 
das Marinhas. Por conta de sua proximidade 
geográfica, a ilha demonstra uma  identidade 
ambiental muito similar a Belém, com 
clima  equatorial superúmido com médias de 
temperatura em torno de 27 "Ce 
pluviosidade de 2.800 mm anuais (Furtado e 
Silva Jr, 2009). Apresenta uma populacáo 
residente estimada emaproximadamente 
27.000 habitantes (IBGE, 2010). 
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Figura 1. Ilha de Mosqueiro destacando em vermelho a localizacáo espacial dos c 
riadouros de caramujo Biomphalaria 


Foi  percorrida uma extensa área para 
demarcacáo e caracterizacáo das colecóOes 
hídricas locais que apresentavam moluscos 
planorbídeos do género Biomphalaria, em 
parceria com os agentes comunitários de saúde 
(ACS) do Programa de Saúde da Família (PSF) 
envolvidos na pesquisa. A  localizacáo 
espacial dos 20 criadouros foi determinada 
através de um receptor GPS (Global 
Positioning System), como pode ser observado 
na Figura 1. As coordenadas dos criadouros 
foram comparadas com os dados dos 
setores  censitários,  obtidos do  IBGE 
(Instituto Brasileiro de Geografía e Estatística). 
Foi utilizado um Sistema de Informacóes 
Geográficas (SIG) para armazenar, organizar e 
analisar as imagens do satélite SPOT 5 e, 
também, para gerar um banco de dados. 

Na construcáo dos mapas de superfície fol 
aplicado Oo kernel nos  criadouros de 
Biomphalaria, segundo a metodologia descrita 
por Guimaraes (2010). A largura de banda de 


300 m foi utilizada para estimar a 
densidade do kernel dos  criadouros. Os 
aglomerados obtidos com o kernel foram 


caracterizados com imagens do satélite SPOT 5. 


3. RESULTADOS E DISCUSSAO 


Todos os bairros pesquisados do DAMOS 
apresentavam  condicóes propícias para 0 
desenvolvimento e proliferacáo dos 
caramujos, pois contavam com presenca de 
corpos hídricos com baixa mobilidade de água, 
abundancia de  vegetacáo  macrofítica e 
contaminacáo de valas e canais por matéria 
orgánica. 

B. straminea foi a única espécie de 
Biomphalaria encontrada nos criadouros após a 
identificacáao morfológica. Porém, alguns 
exemplares seráo submetidos á identificacáo 
por técnicas de biologia molecular (PCR) para 
confirmacáo da espécie. Caramujos de outros 
géneros como Physa e Pomacea, embora 
náo participem do ciclo de vida do S. 
mansonit, foram encontrados com bastante 
frequéncia nas áreas pesquisadas, algumas 
vezes dividindo o mesmo habitat com 
caramujos do género Biomphalaria. 


Tabela 1. Setores censitários, 


bairros trabalhados, 


quantidade de criadouros (PT) e 


porcentagem de domicílios ligados a rede geral de água e esgoto 


Setor 

000001 Vila 

000002 [Vila 

000005  |[Maracajá 
000006  |[Maracajá 
000017 — |[Sáo Francisco 
000022  [Carananduba 
000024  [Caruara 
000035 — |[Sáo Francisco 


Todos os vinte pontos georreferenciados 
apresentaram  colecóes hídricas  contendo 
criadouros. Eles estáo situados em oito setores 
censitários, dois setores pertencem a localidade 
de Maracajá com nove criadouros, dois a Vila 
com trés criadouros, ambos localizados no 
sudoeste da Ilha do Mosqueiro, e na parte oeste 


== IV 





1238 5. O 112,4 
637 89,7 29,7 
1033 950  B1,7 
2312 90,7 5,7 
866 59,7 4,7 
1033 43,6 (0,4 
794 13 0,0 
1572 35,5  [P,8 


da ilha, existem dois setores pertencem a Sáo 


sx 


Francisco com seis criadouros, dois á 
Carananduba com um criadouro e um a 
Caruara com um criadouro. A  Tabela 1 
mostra Os olto setores censitários onde 


foram encontrados criadouros de caramujos 


Biomphalaria. 
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Figura 2. Aglomerados de criadouros de Biomphalaria resultantes da aplicacáo do kernel 


Podemos observar na úTabela 1 que a 
localidade de Maracajá apresenta 92,1% dos 


seus domicílios ligados a rede geral de água 
e 13,7% ligados a rede geral de esgoto. Vila 


tem 93,1% dos seus domicílios ligados a 
geral de esgoto. Sáo Francisco tem 44,4% dos 
seus domicílios ligados a rede geral de água e 
5,4% ligados a rede geral de esgoto. Caruara 
apresenta 1,5% dos seus domicílios ligados a 
rede geral de água e náo tem rede geral de 
esgoto. Carananduba tem 43,6% dos seus 
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rede geral de água e 18,7% ligados a rede 
domicílios ligados a rede geral de água e 
0,4% ligados a rede geral de esgoto. Logo, 
os resultados da Tabela 1 mostram que a 
regiáo oeste apresenta uma maior precariedade 
em relacio ao saneamento, do que a parte 
sudoeste da Ilha de Mosqueiro. 


Figura 3. Aglomerados de criadouros de Biomphalaria com a imagem de 


satélite 


O resultado utilizando o kernel apresentado 
na Figura 2 mostra quatro aglomerados de 
criadouros de Biomphalaria. Um aglomerado 
na regiáo sudoeste e trés na regiáo oeste. 

Na Figura 3 a imagem de satélite mostrou que 
o aglomerado da regio sudoeste (Figura 3A) 
apresenta uma maior quantidade de residéncias 
do que da regiáo oeste (Figuras 3B, C e D), 
principalmente no aglomerado de  Caruara 
(Figura 3B) que apresenta uma área com 
grande possibilidade de expansáo residencial. 


4. CONCLUSOES 


Este estudo mostrou a importáncia da utilizacáo 
de Sistemas de Informacóes Geográficas (SIG) 
e Sensoriamento Remoto para estudar 0 
risco de  esquistossomose na Ilha do 
Mosqueiro, Guimaraes et al. (2010) usou SIG 
para estudar e controlar a esquistossomose em 
Minas Gerais. 

Embora náo tenha sido encontrada qualquer 
Biomphalaria infectado, existe relato de 
moradores com esquistossomose náo autóctone 
na área de estudo. Além disso, os resultados 
mostram que a regio oeste da Ilha de 
Mosqueiro apresenta um alto potencial para 
a Instalagio de focos de transmissáo da 
esquistossomose, pois exibem colecóes hídricas 
com populacáo intensa de caramujos vetores, 
precariedade do sistema de  esgotamento 
sanitário, deficiéncia na distribuicáo de água 
tratada, dentre outros fatores. 

É importante salientar que o DAMOS passa por 
um processo de ocupacáo desordenada de seu 
espaco, com intenso fluxo migratório de 
famílias numerosas que se deslocam para locais 
que náo dispóem de estrutura básica de 
moradia, oriundos principalmente de 
municípios vizimhos como Santa Bárbara e 
Santo António do Tauá, bem como de outras 
regides do país. Destacando os migrantes da 
regido nordeste (área endémica de 
esquistossomose) que vem em busca de 
melhores oportunidades de trabalho, mas que 
em virtude de baixa escolaridade acabam por 
engessar a massa de desempregados, que 
acabam por invadir áreas de matas primárias, 
localizadas em geral próximas a mananciais. 
Por isso, devido ás atividades antrópicas, os 
mananciais tornam-se poluídos, e aliados a isto, 
temos ainda um crescimento populacional 
desordenado de maneira que O sistema de 
saúde local náo tem como atender toda a 
populacáo, contribuindo para O aumento da 
incidéncia de  várias  doencas  infecto 
parasitárias. 

Peixoto e Machado (2005) e  Anaruma 
Filho e Santos (2007) no  estudo da 


esquistossomose na Usina Hidrelétrica de 
Miranda//MG e Campinas/SP mostram os 
mesmos fatores apresentados nesse estudo, 
relacionados aos problemas socioambientais, 
económicos e de falta de servicos essenciais de 
saneamento básico e água tratada. 

A área de risco para a esquistossomose foi 
obtida pelo estimador de kernel que identificou 
os locais através dos aglomerados e, também, 
foi possível calcular, utilizando um SIG, a sua 
extensáo e com o auxilio de uma imagem de 
satélite determinar as características 
ambientais da área do aglomerado. Cardim et 
al. (2011) identificou áreas de risco para a 
esquistossomose em Lauro de  Freitas/BA 
usando dados socioambientais, epidemiológicos 
e estimador kernel. 
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ABSTRACT 


The special fine sediment dynamics on South 
America coastal environment under the 
influence of the Amazon river had been highly 
documented during the last decades. Thus, 
huge  mudbanks migrate along thousand 
kilometers of coast until the Orinoco delta. 
These banks are constituted by an intertidal part 
formed by consolidated mud and with a subtidal 
area of the mudbank composed by fluid mud. 
The aim of this short paper is to combine 
remote sensing and in-situ measurement in the 
determination of mud bank migration both in 
the  intertidal and  subtidal area. The 
determination of the delimitation of the 
intertidal part of the bank are monitored using 
remote sensing images acquired at low water 
level. While the waves attenuation are used as a 
proxy of the subtidal part localization. In fact, 
waves are dampened by fluid mud at proximity 
of the subtidal area, and are highly detectable 
on SPOT images. Moreover,  1n-situ 
measurement in the Kourou river channel 
provide data of the fluidify mud on the channel 
from 2002 to 2012. It was observed that the 
passage of the mudbank contribute differently 
along the channel with a more pronounced 
arrival of fluid mud in the offshore part of the 
navigation channel between 2003 and 2009 
while the onshore channel were silting-up by 
the trailing edge of the mudbank from the mid- 
2009. Combining both the mudbank position 
from remote sensing images analyses and mud 
lens survey allow the determination of the 
localization of a mudbank along the Kourou 
coast for 10 years. 


Keywords: Intertidal, Subtidal, fluid mud, in- 
situ measurement, image Processing 


1. INTRODUCTION 


The Guianese coast 1s the muddiest shoreline of 
the world. These coasts experience important 
fine-grained sediment supply due to the 
Amazon river discharge (Augustinus, 1978). 
These cohesive sediments lead to the 
formation of huge mud banks which migrate 
along the Guianese coast (Froidefond et 
al.,1988;Gardel and Gratiot,2005;Anthony et 
al., 2010) Mudbank and associated 
mangroves-fringe shores form parts of the 
sea side environment along the 1600km from 
the Cabo Cassipore (Brazil) to the Orinoco river 
(Venezuela) (Augustinus et al., 2004, Warne et 
al., 2002). As a consequence of mudbanks 
migration, the shorelines are in perpetual 
fluctuations (Gardel and Gratiot, 2004; Anthony 
et al., 2010). The migration of these huge 
mudbanks (10 to 60km-long, 20 to 30km-wider 
and 5m-thick, Froidefond et al., 1988) could 
attainn up to 5km by year (Gardel and Gratiot, 
2005). This coast is organized in “bank” and 
in “interbank” phases(Anthony and 
Dolique,2004; Baghdadi et al.,2004). Bank 
phases consist of a prograded coast protected by 
wave action due to the important wave 
dampening on the mudbank offshore extension 
(Anthony et Dolique, 2004). 


Bank phases lead to the installation of 
mudbank, and therefore a prograded coast from 
20 to 30km shore to offshore, consolidation of 
mud  engendered the development of a 
mangrove forest. Besides, these phases lead 
to coastal silting-up no matter the type of 
environment (sandy beach, rocky shore, harbor 
or harbor channel). As a contrary during 
“Interbank” phases the coast is free of any 
protection against wave attack, and erosion of 
the colonized mudbank and shoreline are 
consequent. Mudbank are divided in several 
part, the first one consist of an intertidal part 
which 1s situated in the fluctuated tidal 
environment. This part 1s composed by 
consolidated mud created by fluid mud arrival 
through the action of coastal current. The 
subtidal part which is in permanent submerge 
are formed by fluidify mud, the waves are 
totally attenuated at their approach. The subtidal 
area of a mudbank constitutes the bigger part of 
the bank and therefore their migration represent 
a key in the migration processes. Several 
studies had been focused on intertidal mudbank 
migration (Froidefond et al., 1988, Gardel and 
Gratiot, 2004,2005) with the purpose of 
quantified the migration rates of a mudbank 
along  interbank area. Nevertheless, any 
geomorphological studies had been devoted on 
the characterization of the subtidal part 


migration processes, which need important in- 
situ measurement such as seismic profiles and 
bathymetric survey. 

The extreme coastal fluctuation along the 
French Guiana coast engendered important 
difficulties for coastal decision makers. A 
better knowledge of the integral mudbank 
migration processes 1s then a priority in order to 
have a better understanding of the rapid coastal 
changes. The aim of this paper 1s to integrate 
both data extracted from remote sensing 
analysis and in-situ measurement for tracking 
mud-bank migration. 


2. METHODOLOGY 


2.1 Study site 


The study site is located on the Kourou coast, 
which comprised the Kourou river channel 
harbor (e.g .Figure 1). These coast suffer from 
important fluctuations due to the complex 
mudbank processes. Arrival of fluid mud on 
the coast leads to the modification of coastal 
bathymetric. 





Figure l. Satellite imagery (SPOT image, 30/09/2011, SEAS Guyane) of the Kourou coast and 
localization of the Kourou river 


Along Kourou coast, arrival of mud cause a 
silting-up of the coastal environment and could 
lead to a silting-up of the Kourou river 
navigational channel, provoking the non-access 
of the river due to the high elevation of the 
bottom depth. Therefore, constant dredging are 
realized to evacuate the mud inside the channel, 
and to facilitate access of the navigational 
channel. 


2.2. Remote sensing 
To determine coastal change, multi-spectral 


and multi-temporal satellite images are used in 
the present study. Intertidal part of mudbank is 
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clearly visible on satellite images when 
acquisition is realized at lower water level, 
whereas subtidal part could be delimited only 
by waves dampening. This type of data permits 
a geo-localization of mudbank on the Guianese 
coast. Satellite images are used to extract the 
position and area of the mudbank using a semi- 
automatic algorithm (e.g. Figure 2.A). This 
method  regroups  photo-interpretation and 
segmentation. The two-dimensional 
morphology of the mudbank is obtained at 
lower water- level. 
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Figure 2. Extraction methodology: A/two-dimensional and B/three-dimensional morphology 


The near infra-red (NIR), mid  infra-red 
(MIR), red and green bands are used to 
extract the coastal features elements of our 
study: mangroves-fringe and mud. 

The surface of intertidal part of mudbank 
depends on the water level as a consequence 
images acquired at different water-level allow 
the realization of three-dimensional mudbank 
morphology. 

Therefore, the intertidal topography of a 
mud-bank is obtained by several isocontours 
acquired at different water level between low 
and high tide (e.g. Figure 3). 

Two conditions are required for quality results: 
images acquired in a short-time interval (to 
analyze the same bank position) and a panel of 


images which cover the tidal amplitude (e.g. 
Figure 3). These data lead to the realization of a 
Digital Elevation Model (DEM) of the muddy 
bank. 

The interval of altitude of the mud correspond 
to the interval between two isocontours. All 
images were in the world reference 
system(WGS 84) datum and were projected 
using the Universal Transverse Mercator system 
(zone UTM 22). 

The subtidal part of the mudbank had been 
extracted by manual method on eight SPOT 


images in this study (14/10/2001, 
02/11/2002, 15/09/2003, 25/06/2006, 
24/09/2007, 02/07/2008, 16/10/2010, 
26/10/2011). 


Unte | Water Level | Satelite 4 Resolution | 
10/09/2010 |  0m89 |  SPOTS-10m | 
OS 09/2010 | 17 1 SPOT 5 10 m 
16/10/2010 |  2m82 SPOTS -10 m 
01/10/2010 | 2m40 SPOT 5 -10 m 





Figure 3. Chosen SPOT-5 images for the three-dimensional mudbank morphology 
(Example of the mudbank of Kourou in 2010) 


2.3. In situ measurement 

From 2002, Atlantic Dragage society has been 
mandated by the Guianese Spacial Center with 
the purpose of prevents silting up of the Kourou 
river by monitoring the fluctuated channel 
depth. The data acquired during the bathymetric 
survey permit extraction of the mud lens 
thickness along the Kourou river channel (14- 
kilometers long) on particular point, named 
“Kilometric Point (PK). The mud lens 
thickness 1s an indicator of the potential 
arrival of fresh fluid mud on the channel. The 
long-term analyze of the bottom thickness of 
the mud- lens gives temporal evolution and 
migration of the mudbank. 


3. RESULTS AND DISCUSSION 


3.1. Tracking mud-bank from remote sensing 
Several methods are used for coastal 
environments monitoring of French Guiana. 
However, for a better knowledge of the coastal 
processes involved in the mudbank migration, 
we need the combination of remote sensing and 
In-situ measurement. The present study 
combined several approaches in order to 
evaluate the migration of a mudbank along the 
French Guiana coast. The migration of the huge 
mudbank along the coast modify the shore 
morphology through time. 


The knowledge of their two and  three- 
dimensional morphology and processes of 
migration are required to prevent or to predict 
and follow their movement. 

In 2010, the two-dimensional morphology of 
the mudbank of Kourou (e.g. Figure 4A) 1s 
obtained between coastline (green line) and an 
isocontours acquired at a water level of 0.89m 
(red line). This low level had permitted to 
observe a huge surface of intertidal mudbank. 
A diachronic survey will permit to obtain two 
and three-dimensional deformations of the 
intertidal part of the mudbank using spatial 
indicator (e.g. Figure 4A and 4B). 

In 2010, the DEM of the mudbank of 
Kourou(cf.fig.4) indicated lower altitude on the 
offshore part of the mudbank whereas the 
inner part is more concentrated in mud and 
therefore has higher altitude. Besides, we noted 
that the trailing and lee edge of the mudbank 
comprised lower altitude and therefore less mud 
than the middle part of the mudbank which 
attains up to 2m40. 

The two and three-dimensional morphologies 
are used in order to define deformations: 
migration velocity (distance between two 
successive barycentre, AVB, or ARB), and 
altimetric deformations (difference between two 
DEM). 





Figure 4: A/ Two-dimensional morphology defined by spatial indicators (front of mudbank (AVB), back 

of mudbank (ARB) and barycentre of mudbank). The intertidal mudbank is also characterized by and 

area which is defined between an isocontours acquired at low tide and coastline. B/ Three-dimensional 

morphology; result of a spatial interpolation, two main types of interpolation (kriging and triangulation 
method). (Images: SPOT-3 images from SEAS Guyane) 


The analyze of the subtidal part indicated a 
northwestward migration and a deformation of 
the bank along the coast (e.g. Figure 5). The 
offshore limit of the bank is curved from the 
lee to the trailing edge. The recent subtidal 
delimitation (from 2006) after the passage of the 
Kourou river suggested a more larger subtidal 
part, facilitated by the arrival of recycling mud 
from the back to the front of the bank. This 
observation coincides with the previous one 
for the intertidal part of the mudbank. The 
subtidal delimitation of the mudbank gives the 
total cover area of the bank along the coast. If 
the bank position identification was realized 


only by the intertidal bank there will be a lack 
of data. In fact, by observing the subtidal 
delimitation in 2010 (e.g. Figure 5) and the 
intertidal area (e.g. Figure 4) we could denote 
that the intertidal bank represent a thinner area 
of the total bank. Contrary to the front of the 
intertidal area (e.g. Figure 4) which represents 
the most thinner part of the intertidal zone, the 
front of the subtidal part is the more larger area 
of the mudbank. Therefore, the front of the 
mudbank constitute a bigger source of fluid mud 
during the migration processes. 
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Figure 5. Localization of the subtidal area of the mudbank along the Kourou coast from 2001 to 2011, 
localization of the navigational channel of the Kourou river PK 5800 (at the mouth of the Kourou river) to 
PK 13800 (end of the channel) 


3.2. What happens in the Kourou river? 
Regular bathymetric surveys of the Kourou 
harbor channel offer a detailed situation of fluid 
mud arrival along the channel and so on the 
subtidal migration of the mudbank. Fresh mud 
arrival in the Kourou channel  induces 
significant depth addition, more precisely the 
mud lens, furthermore important dredging effort 
in the navigation channel is organized to allow 
the passage of economical shipping. This mud 
lens represents a significant layer of fluid 
mud which overlay the consolidated bed. In fact, 
fluid mud is easily transported through the 
coastal current and could be deposited on other 
layer of consolidated mud. The bathymetric 
survey highlights the displacement of the bank 
along the navigation channel during the 
measurement period. 

During the study period, sporadic events of 
important arrival of fresh mud in the Kourou 
harbor channel were observed (e.g. Figure 6). 
The lack of hydrodynamic data near this coastal 
site do not permit a coupling between the 
sporadic silting-up and eventual changes of 
the Hhydrodynamic conditions, and thus 
explanation of these arrival of fluid mud only 
could be coupling by the position of the 
mudbank. Heterogeneity of the mud thickness in 
the channel is denoted depending on the 
localization across the harbor channel, more 


precisely the distance from the coast. 

Silting-up of the channel was first observed 
from the mid-2004, the dredging effort were 
then concentrated from the PK 8000 (approx. 
mid-channel) to offshore, in order to liberate 
the navigational channel. Generally between 
June 2002 and April 2010 (except in the 
beginning of 2008), thinner mud thickness 
were observed at the mouth of the Kourou 
river (near PK 5800). 

Between June 2002 and July 2004, more 
important mud  thickness were observed 
between PK 7800 and PK 9800 oscillated from 
0.73m to 2m. The extremity of the 
navigational cannel encountered less arrival of 
fluid mud as indicated by the in-situ 
measurement. Between November 2007 to June 
2008, important fresh arrival of mud were 
observed in the internal channel from PK 5800. 
From 2001 to 2004, the curving subtidal part, 
and therefore the more important part of the 
bank, was localized in front of the Kourou 
channel, and  relatively less mud was 
encountered in the channel(e.g. Figures 546). 
Between 2003 and 2006, the mudbank was 
crossing the Kourou river, the major part of the 
bank was beyond the channel, we noted a 
stagnation and a slowing down of the bank 
migration, as the Kourou river constitute a 
natural barrier. As a contrary, since 2006 


when the major part of the bank were 
situated after the Kourou river channel, it 
appears more faster migration of the mudbank. 

The beginning of the channel river (from PK 
5800 to 6500) seemed to be silting-up since 


2010 (e.g. Figure 6) while the trailing edge of 
the Kourou bank was localized at the vicinity of 
the cannel (e.g. Figure 5). In fact, more 
consistent mud were found in this part of the 
channel during this period (2010 to 2012). 
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Figure 6. Mud lens thickness in the Kourou river channel between 2002 and 2012. PK 5800 indicate the 
position at the mouth of the Kourou river, the navigation channel end at the PK 14000 
(not represented in this figure) 


4. CONCLUSIONS 


Mudbanks influence on the coast require intense 
monitoring in order to evaluate their spatial 
localization on the coast either during the bank 
or inter-bank phases. The study of the nearshore 
processes of mudbank had to be realize at 
different spatial and time-scale. 

The first tools contribute to the localization of 
the bank with two and  three-dimensional 
morphology characteristics such as an idea of 
the altitude of the mud, and the quantity of 
muddy sediment on the mudbank. While in-situ 
measurement is a support to remote sensing 
extracted data by consolidate the results 
obtained by satellite  imagery by the 
quantification of the amount of mud present 
on the shore or estuary at the different stages 
of the mudbank migration. 
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... y Envisat sigue observando el lugar 
del próximo simposio (arriba izquierda) 
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